On a generalization of the Ito-Wentzell formula for 

system of generalized Ito's SDE and the stochastic 

(N ; first integral 

o 

D 

Q I Elena V. Karachanskaya 

^ ■ Pacific National University, Russia, Khabarovsk 

December 6, 2011 

Oh 

' Abstract 
> . 

^ ' Generalization of the Ito-WentzcU formula for the generalized Ito's 

SDE (Ito's GSDE) system with a non-centered measure is constructed 
on the basis of the stochastic kernel of integral transformation. The 
Ito's GSDE system for the kernel the solution of which is the kernel of 
the integral invariant, is formed. This invariant is connected with the 
solution of the Ito's GSDE system non-centered measure. The concept 
of a stochastic first integral of the Ito's GSDE system with non-centered 
measure is introduced and conditions that when being performed the 
random function is the first integral of the set Ito's GSDE system are 
defined. 

Introduction 

^ ' The specific approach to the stochastic integration based on the fact that not 



(N 

O 
(N 



the point but its neighborhoods is reflected according to some quite different 
from the traditional rule, based on Malliavin's calculus, was presented in 
[1, 2, 3]. 

The transfer to the stochastic equations for definite realizations in based 
on comparing the dynamics of the space structures element marked by an 
integer to the dynamics of the trajectory's realization for which this integer 
is the initial value. This possibility is stipulated by the proximity of the 
trajectories, in the sense of probability, beginning the domain of the initial 
neighborhood for each of the temporary intervals. 

The investigation of open systems that are characterized by the presence 
of random disturbances in them both continuous described in Wiener's pro- 
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cesses and jump-like ones, presented in Poisson's processes is urgent nowa- 
days. At the same time in case of random disturbances some characteristics 
of the system must be presented. And the problems of constructing the pro- 
gram control in the open systems arise. When a disturbance is caused by 
Wiener processes only there are some types of algorythms for constructing 
program control. The Poisson component of the random disturbance still 
not let to do program control for such systems. The present generalization 
of the Ito-Wentzell formula solves this problem [5]. 



1 The stochastic volume 



Let x(t) ~ is the random dynamic process and this process is solution of the 
next SDE system: 

m ~ 

dxi(t) = ai (^t; x(t)^ dt + '^ bi^k (t; x(t)^ dwk(t) + gi (t; x(t); 7^ ^7), 



k=l 

x(t) =x(t;x(0))|^^^ = x(0), 



i = 1, n, 



t > 0, 



where w(t) is m~dimensional Wiener process, A7) - is homogeneous on 

t non centered Poisson measure. 

We suppose that the coefficients a{t; x), b{t; x), and g{t; x; 7) are bounded 

, 111 9aj(t;x) dbik{t;:x.) dgAt^n.;-^) 

and continuous together with own derivations 



dxj 



dxj 



dxj 



i,j = l,n respect to (t;x;7). 

These conditions are sufficient for existing random variables 



Jij{t) = 1 .i.m. 



Xi(t;x(0) + A) -Xi(t;x(0) 



that can be defined by solving the system (1) and the system of stochastic 
equation [4]: 



dJijit) = Jij{t) 

+ 



dxi{o) 

dg^{t-^{t)) 
dxm 



v{dt\ d^) 



t=0 



0, i^j, 

1, i=3- 
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Let J(t) = J(^t;x(0)j = det Jj^ (^t; x(0)j j , i,j = l,n. Then the ran- 
dom variable J{t) can be considerer a Jacobian of transition from x(0) to 
x(t;x(0)). 

By means of random variables Jij{t) it is possible to construct different 
random manifolds M (t; j^{t; X)J (A G P C W, x(t; A) G M", r < n), 
having compared them to configuration manifolds of the initial manifold 
M f 0; x(0; A)) (X € V C R\ x(0; A) = x(A) G M", r < n). 



Definition 1 [2] Let's call the structure dr{t) = J{t) dT\^x{0)j , where (ir(x) 

n 

JJdxj, an element of stochastic volume that it generated by some random 

i=l 

process x^t; x(0; A)^ , where J(t) is a Jacobian built of elements Jij (^t; x(0) 
According to definition 1 a random structure is introduced: 



Xr(t)(i) = f{t;z)dT{z) = j ■■■ j /(t; x(t; y)) dr(y), 



(3) 



T(t) r(o) 



where x(i; y) is a stochastic process with random initial condition y. 

in so doing the integration over a stochastic volume T{t) is carried out 
Riemann's integration sums, mean square limit concept for any continuous 
functions by z, and in general, random functions /(t;z). 

Generically, the structure (3) is a random variable. The statement of 
the classic integral calculus can be associated with random structure (3). In 
this sense the equality (3) is an analogue for substituting variables in the 
classic calculus. Here, a Jacobian of transformation Jij{t) is a random value 
comparising the solutions of GSDE system (2). 

Let us determine the conditions for the invariance (or conservation) of 
the stochastic volume. 



2 Stochastic kernel of the stochastic integral in- 
variant 

Let's consider the Ito's GSDE system of the following kind (1). Let's assume 
that p{t; x; tj) - is a random function which is measured with respect to cr- 

algebras flow , J'^ti C J-'t2, ti < t2. This flow is adapted with the 



3 



processes w(i) and i^(t;A7). Each function /(t;x) from the set of local 
bounded functions which have second bounded derivatives with respect to 
X the next relations hold: 

J p{t; x)/(t; x)dr(x) = J p(0; y)/(i; x(t; y))dT{y) (4) 
p(0;x)dr(x) = 1, (5) 



lim p(0;x) = 0, dT{x) = Y\dxi. (6) 

x|— )>oo 

1=1 

Here x(t;y) is the solution of (1). 

If we have f{t;x) = 1, then taking into consideration the conditions (4) 
and (5) we'll have the equality 

I p(t;x)dr(x) = I p{0;y)dT{y) = 1. (7) 

M" R" 

It means that there exists a random functional, which retains the constant 
value for the random function p{t;x) = p{t;x;u}): 

J p{t;x)dr{x) = l. (8) 

We consider this functional a stochastic volume. Then the equation (4) 
under conditions (5) and (6) is a stochastic integral invariant for function 

/(i;x). 

Definition 2 Let's call the non-negative function p{t;x) a stochastic kernel 
or stochastic density of the stochastic integral invariant if equation (4), (5) 
and (6) are kept. 

Remark 1 But the essential distinction that let to consider the invariance 
of the random volume on the basis of integral invariant kernel in [2, 3] to 
that there is a functional factor in (4) . Thus, the notion of integral invariant 
kernel for the system of the ordinary differential equations can be considered 
a particular case of the introduced one in [2, 3] if we take /(i;x) = 1 and 
examine the integration by determined volume, having excluded the random- 
ness preset in Wiener and Poisson processes. 



4 



Let arbitrary function /(t;x) has the first derivation with respect to t 
and the second derivations wish respect to x. We determine relations for 
function p{t;x) which imphes that p{t;x) is integral invariant kernel for 

/(t;x). 

Let x(t) - is solution of (1) and /(t;x(t)) is random function. We apply 
the generalized Ito's formula [4]: 



<i,/(t;x(t)) 



i=l 



+ 



n n m 



i=i j=i k=i 

n m 



dt+ 



+ 



i=l k=l ^* 

[/ (t; x(t) + g{t; x(t); 7)) - /(t; x(t))] u{dt; d-f). 



(9) 



Remark 2 Later on for notational simplicity we 'II not use the summation 
symbol if we have double occurred indexes. 

Let's differentiate both parts of equation (4) wish respect to t taking 
into account that /(t; x) of the left part is a determined function while 
/(t;x(t;y)) of the right part is a random process. 

/(t;x)dip(t;x)+Ki;x)^^^^dt)dr(x)= / p(0; y)(ii/(t; x(t; y))dr(y) 



Then, taking into account both (4) and (9) we get: 



f{t; x)dtp{t- x) + p{t- x)^i^dt ) dr(x) 



p(0;y)dj(t;x(t;y))dr(y) = / /^(t; x)dt/(t; x)dr(x) 
dr(x)/5(t;x) 



9/(t;x(t)) , 9/(t;x) , 
+ ai(t;x) — h 



9t 



1, / N, / ,d^f(t;x) 
+-6.,(t;x)6,,(t;x)^-^ 



+ 



+ 6i fc (t ; x) dw;fc (t ) + 



/(t;x + g(t;x;7)) -/(t;x) u{dt;d-f) 



(10) 
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Let's consider the integral 



/ = y (ir(x)p(t;x)/(t;x + 5(t;x;7)) 
We change variables: 

x + 5'(t;x;7) = y, 



fill 



(12) 



Having taken into account the transition from x to y, we get x = y 
g(t;x]j) = y — ^(t; x~^(i; y; 7); 7). And the integral / becomes: 



1 = 1 dr{y)p {t; y - g{t; ^-\t; y; 7); 7)) • f{t; y) • D {^~\t; y; 7)) , 



where D (x ^{t; y; 7)) is a Jacobian of transition from x to y. 
Taking into accout (13) and integration over the we have: 



(13) 



dr(x)p(t;x) J ^/(t;x + 5(t;x;7)) -/(t;x)Ji/(dt;d7) 

IR(7) 

dr(x)p(i;x) j f{t-x + g{t-x;^))u{dt;dj)- 

R(7) 

c?r(x)p(t;x) j f{t;x)u{dt;dj) = 

18(7) 

y"dr(x) j (p(t;x-<7(t;x-i(t;x;7);7)) • 



(14) 



•/(t;x).L»(x-i(t;x;7)) jz.((it;(i7)- 
dr(x)p(t;x) y f{t;x)uidt;dj) 



dr(x)/(t;x) y [p(t;x-5(t;x-i(t;x;7);7)) 

R(7) 

•Z)(x-i(*;x;7)) -P(i;x)li/(dt;d7)- 



With due regard for (5) we calculate the following integrals using integration 
by parts: 

5/(t;x) 



dT{x)p{t;x)ai{t;x)- 



dxi 



(15) 
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dT{x)p{t;x.)bikit;x 



5/(t;x) 



dT{x)p{t; :K.)bi k{t; x.)bjk{t; x) 



dxi 
52/(i;x) 



dxidxj 



(16) 
(17) 



For equation (15) we get: 

dT{x)p{t;:>c)ai{t;x.)^^^^ = f dxjp{t; x)ai (t; x) ^"^^^ "^^ 



+ 00 



dxj J p(t;x)oi(t;x) 

J=i -oo 



g/(t;x) 

dxi 



dxj 



Then we calculate the inner integral using integration in parts, having taken 
into account (6): 



+ 00 



p(t\yi)ai(t\ x) — dxi 



dxi 



p{t;x)ai{t;x)f{t;x) 



+ 00 



+ 00 



d{p{t;x)ai{t;x)) 
dxi 



dxi 



/(t;x)^M^i^Mi^d^,. 
dxi 



Therefore, we get 



dT{x)p{t;x)ai{t;x) 



df{t;x) 

dxi 



+ CXD 



dxjp{t; x)ai {t; x 



i=i 



dfit;x) 
dxi 



+ 00, 



'n-1 



+ 00 



-oo 

+ 00 



n 



Y\dxj / p{t;x)ai{t;x) 



dfit;x) 
dxi 



dxj 



-oo 
-l-oo 



\{dxj / /(t;x) 

i=i -oo 



(ir(x)/(t;x; 



d{p{t;x)ai{t;x)) 
dxi 

d {p{t]x)ai{t;x)) 
dxi 



dxj 



(18) 
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Similarly we calculate the second integral (16) and the third integral (17) 
by using twice the integration in parts: 



dr(x)p(t;x)6ifc(t;x) 



g/(t;x) 



dp{t;x)bikit;x.) 



dxi 



(ir(x)/)(t; x)6j fc(i; x)6j ^ {t; x) 



gV(t;x) 

dxidxj 



5 2 {p{t;x)bikit;x)bjkit;x)) 



(19) 
(20) 



dxidxj 



In equation (10) we move all terms to the right part and taking into 
account (14), (18), (19) and (20) we obtain the following equation: 



= j dr(x)/(t;x) 



9/(t;x(t)) , 
-dt (/o(t;x)) + 



a/(t;x(t)) dp{t;x)bik{t;x) 



dt 



dxi 



dt 



dwk{t)+ 



+ 



d{p{t;x)a,{t;x{t))) , 1 5^ (p(t; x)6ifc(t; x)fej fc(t; x)) 



+ 



+ 



dxidxj 



dt+ 



dxi ' 2 

P - fif(t;x"^(t;x;7);7)) • L» (x"^(i; x; 7)) -p(i;x) v{dt;d'y) 



(21) 

To make the equation (21) true for any locally bounded function f{t;x.) 
obtaining the limited second-order derivatives the integral invariant p{t; x) 
must be the solution for the following stochastic equation: 

d,pit;x) = + ^_ 



dxi 

^ I {p{t;x)bik{t;x)bj kit;x.)) 
2 dxidxj 



dxi 



dt+ 



+ J [^(t;x-g(t;x-i(t;x;7);7)) •Z)(x-i(t;x;7)) -p(t;x)]i/(dt;d7). 

18(7) 

(22) 



The following terms must be kept: 



p{t;x) ^^^=p(0;x) e Cl 
lim p(0; x) = 0, lim = 0. 

x|— >oo |x|— >oo OXi 



(23) 
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Thus, we get the invariance conditions of a stochastic volume and the fol- 
lowing theorem is proved. 



Theorem 1 Let x{t), x G R" - is the solution of ltd' s GSDE system: 

m „ 

dxi{t) = ai (^t; x{tfj dt + '^ bi^k (t; x{t)^ dwk{t) + j Qi (t; x(t); 7^ i^{dt; d'j), 

M(7) 

x(t) = x(^t;x(0)^ 1^ ^ = x(0), i=T~n, t > 0, 

where w(t) - is m-dimentional Wiener process, v{t] A7) - is homogeneous 
on t non centered Poisson measure. Assume that p{t; x) is a random func- 
tion which is measurable function with respect to a-algebras flow |-^| ; 

Ttj^ C J-t2, ti < t2, which adapting w(t) and A7) processes and with 
respect to any function f{t;x) € © from the set of locally bounded func- 
tions having the second bounded derivatives in x. The function p{t; x) is the 
stochastic kernel of the stochastic integral invariant (4) for arbitrary locally 
bounded function /(t;x) € S if it's a solution of Ito's GSDE system (22) 
meeting the initial conditions (23). 



3 The generalization of Ito-Wentzell formula 

The equation (4) under conditions (5), (6) and equation (22) for the integral 
invariant kernel conformed to a determinate function /(t;x). Let's assume, 
that analogical equation for random function z(t; x) holds: 

j pit;x)z{t;x)dT{x) = j /)(0; y)z(t; x(t; y))dr(y), ^24) 

here x(t; y) is a solution of (1). 

Let's consider a compound stochastic process z(t;x(t;y)) G where 
x(^i y) ^ is a solution of (1) where the process z(t;x) - is a solution of the 
system: 



dt2{t;x)=Yl{t;x)dt + Dk{t;x)dwk{t)+ j i^{dt; d'y)G{t;x;-f). 



(25) 



As far as the random functions n(t;x), Dk{t;x), G{t;x;j) defined on 
are concerned we can assume that they are continuous and bounded 
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together with their own derivatives with respect to aU variables. These 
functions are measurable with respect to cr-algebras flow J-^, which adapted 
with processes w(t) and z^(t; A7) from equation (1). 

Relying on the equation for the stochastic integral invariant. We can con- 
struct the differentiation law for the compound stochastic process z (t; x(i; y)) . 

Let us consider the integral J p{t;x)z(t;x)dT{x.). As we do integration 

over the M" and x € M"", we use equation (24) and having changed the places 
of the equation, we get : 

J p(0;y)z(t;x(t;y))dr(y) = J x)z(t; x)dr(x). ^26) 

Lrt's differentiate both parts of (26) with respect to t: 

p(0;y)(itz(i;x(t;y))dr(y) = 

I (p(t; x)dtz{t; x) + z(i; x)dtp{t; x) - D,{t; rft) dT{x), 

R" 

(27) 

Taking into account that integration is carried out over the and the 
variable is changed we should keep in ind this substitution (12) when inte- 
gration over the curve R{'y)- 

As the function p{t; x) - is the integral invariant kernel (24) we apply 
the theorem 1. Substitute (22) and (25) to the right part of (27): 

h = f (p(t; x)dtz{t- x)+z(t; x)dtp{t; x)-ZJfc (t; x) "'j^^/' dt) dr(x) = 



dr{x) I p{t; x) \U{t; x)dt + {t; x)dwk (t) + 
+ j G{t;x + g{t;x-\t;x;-f)-j)u{dt;dj) 



+ 



+z{t;x) 



dp{t;x)bik{t;x) 
^- dw,{t)+ 



d{p{t;x)ai{t;x)) ^ 1 9^ {p{t; x)bik{t; x)bjk{t; x)) ' ^ 
V dxi 2 dxidxj 
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+ 



p(t;x-c/(i;x ^(t;x;7);7)) • I? (x ^(t; x; 7)) - x) u{dt;d-f) 



-D,(t;x)M!liSM!liS.rfi^. 



Collect similar terms: 



Ii = J dr(x)(^p(t;x)n(t;x) -z(t;x) 



d {p{t;x)ai{t;x)) 
dxi 



+ 



-Dfc(t;x; 



'9p(t;x)6^fc(t;x) 1 (p(t;x)bifc(t;x)bjfc(^;x)) 



5Xi 



dxidxj 



dt+ 



+ I dr(x) x) x) - z(t; x) ^ ^j^^^^^^^' ) dw,it) + 
+ j dT{x) p(t;x)- y G(t;x + (7(t;x-i(t;x;7);7))z.(di;d7) + 

K" /?(7) 

+z(t;x) y p (t;x - c/(t;x~^(t;x;7);7)) • D (x"^(t;x;7)) - p(t;x) i/(dt;d7) 

(28) 

And because of (18), (19), (20) we have: 



(ir(x)z(t; X 
dr(x)Z)fc(t;x 



g (/o(t;x)aj(t;x)) 
dxi 

d{p{t;x)bik{t;x)) 



dxi 



dT{x)p{t;x.)ai{t;x) 
dr{x)p{t;x)bik{t;x 



dz{t; x) 



dxi 

dDk{t;x) 
dxi 



j c?r(x)z(t;x 



{p{t;x)bik{t;x)bjkit;x)) 



dxidxj 



dr{x)p{t; x)bi k {t; x)bj ^(i; x) 



9^z(t;x) 
dxidxj 



Calculate the last integral in (28): 

h = J dT{x)z{t;x) ■ p(^t;x- g{t;x~'^{t;x;-f);-f)^D{x~'^{t;x;-f)). 1^29) 
ffi" 

Let's introduce the variables' substitution: 

x-g(t;x-^(t;x;7);7) = y, 
X = y + g{t; x~'^{t; x; 7); 7) = y + g{t; y; 7). 
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Let L)o(x^^(t; y; 7)) be a Jacobian of transformation from the vector x to 
the vector y. Then, by virtue of (11) and further formal substitution of the 
integration variable's notation, we get: 

l2= j dr(y)z(t;y + 5(t;y;7)) ■p{t;y)Do{x-\t-y--f))D{^-\t-^-j)) = 

= J (ir(y)z(t;y+(7(t;y;7))p(t;y) = J dr(x)z(t; x + ^(t; x; 7)) x). 
MP- R" 

(30) 

As a result, the right part of (27) becomes: 



Ii= J dr(x)p(t;x)|(z)fc(t;x)+5ifc(t;x) 



dz{t; x) 

dxi 



dwk{t)+ 



+ (n(t;x) + a,(t;x) ^' ^ +6,fc(t;x) t} ' ' + 

\ OXi OXi 

1, / N, / xC^^z(t;x)\ , 
+-6..(t;x)6,,(t;x)^-^)dt+ 

+ y" G(t;x + 5(t;x-^(t;x;7);7)i/((it;(i7) + 



(31) 



R{i) 



z(t;x + 5r(t;x;7)) -z(t;x) v{dt]d'y) 



+ 
«(7) 

In equation (27) we transfer all the terms to the right part. And having 
taken (4) into account, we have: 

= ydr(y)p(0;y)|-(itz(t;x(t;y)) + 

+ [Dk{t; x(t; y)) + bik{t] x(t; y))M^^ dwk{t) + 

/t-t/ / xn / x5z(t:x) , , .dDk(t::x.) 
+ n(t; x(t; y)) + ai(t; x) ^ '/ + ^^^(t; x) ' ' + 

+ h,,{t; x(t; y))6, ,(t; x(t; y))^Mg|iZll) dt+ 

+ j G(t;x(t;y) +9(t;x(t;y);7)i^(dt;d7) + 

+ j [z(t;x(t;y) +5(t;x(t;y);7)) -z(t;x(t;y))]zy((it;(i7)|. 

i?(7) 
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Therefore, we'll have the differentiation law for the compound stochastic 
process as the theorem. 



Theorem 2 Let x(t; y) G M" is the solution of SDE system (1), z(t;x) 
is the solution of SDE system (25) and z(t;x(t;y)) € R"° is the random 
process. The random functions that are coefficients of equations (1) and 
(25), which are defined on and M" respectively. This functions are 
continuous and hounded together with their own derivatives respect to all 
variables. Also these functions are measurable respect to a-algebras flow Tt, 
which is adapted with processes w(t) and iy{t; A7) from equation (1). Then 
the compound stochastic process z(t;x(t;y)) is the solution of SDE system: 

dtz{t;^{t;y)) = (Dk{t; x(t; y)) + 6,fc(t; x(t; y)) ^^^^'gj^' ^^^ )dwk{t) + 



n(t;x(t;y))+a,(t;x(t;y)) 



1 



az(t;x(t;y)) 
dxi 



+-bik{t; x(t; y))bjk{t; x(t; y)) 



+ 



+ ^jfc(i;x(t;y)) 
52z(i;x(t;y)) 



dDk{t;x{t;y)) 



dxi 



+ 



dxidxj 



dt+ 



+ 



G{t- x(t; y) + g{t- x(t; y)-^)v{dt; d^) + 
z(t;x(t;y) +9(t;x(t;y);7)) -z(t;x(t;y)) u{dt;d-i). 



(32) 



On the analogy of the well-known terminology let's call the formula 
(32) as generalized Ito-Wentzell formula for the generalized SDE Ito of non- 
centered Poisson measure. In addition we can settle one more statement. 



Proposition 1 // z(t; x(t; y)) - is the solution for equation (32) meeting 
the initial condition 

z(i;x(t;y))[^^ = z(0;y), z(0;y) € Cl 

then the random function p(i;x), which is a function that measurable with 

respect to a-algebras flow |-^| > -^h C /t2> *i < ^2, that is adapted with 
processes w(i) and iy{t; A7), is the stochastic kernel of the stochastic integral 
invariant (24). 
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4 The stochastic first integral 



In the article [1] the concept of the first integral for Ito's SDE system (with- 
out Poisson part) has been introduced. In the article [3] the concept of the 
stochastic first integral for Ito's GSDE system with centered measure has 
been given. Let's introduce the similar concept for the non-centered Poisson 
measure. 

Definition 3 Let random function u{t;x;uj) and solution x(t) of the Ito's 
GSDE system (1) are defined on the same probability space. The function 
u{t; x; ui) is called the stochastic first integral for Ito's GSDE system (1) with 
non centered Poissin measure, if condition 

u(t;x{t;yL{0));u?j = u(0;x{0);u?j 

is held with probability is equaled to 1 for each solution x(t;x(0);a;) of the 
system (1). 

We set conditions for the function u{t;x;uj) which imply that u{t;x]uj) 
is the stochastic first integral for Ito's GSDE system (1). 

Lemma 1 // the function p{t; x) is the stochastic kernel of the integral in- 
variant of order n of the stochastic process x(t) outgoing from a random point 
x(0) = y, then for any t it complies with the equation p^t; x(i; y)^ y) = 

p(0;y), where J(t;y) = J(t;x(0)) is a transition Jacobian from x(t) to 
x(0) = y. 

Proof. Let's proceed to the equation (7). If we substitute the variables 
of the integral we'll get the integration is fulfilled according one and the 
same random volume: 

1 = I p{0;y)dr{y) = J p(t;x)dr(x) = J x(t; y)) J(t; y)dr(y). 

H" M" M" 

As the random process x(t) is defined within the extended phase space of 
general dimension M""*"^, then to determine the uniqueness of the trajectory 
the system of differential equations for stochastic kernel should consists of 
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at least n + 1 equations: 

OXi \ OXi 

^ I d^ {pi{t;x)bii,{t;x)bjkit;x)) \^ ^ 
2 dxidxj J 



pi{t;x{t)) = pKO;x(0)) = pz(0;y), l = l,n + l. 

(33) 

As you know the collection of kernels is called complete if any other 
function that is the integral invariant kernel of the n-order can be presented 
as the function of this collection' elements. 

Theorem 3 The GSDE system (1) possesses the complete collection of ker- 
nels, which consists of {n + 1) functions and each of these functions is the 
solution for equation (22). 



Proof. Let p;(t;x), / = l,m, m > n + 1 - are the integral invariant 
kernels of (4). As follows from lemma 1 that for any I ^ n + 1 the ratio 

— ^ ' \ is a constant depending only on x(0) = y for any solution 
Pn+i(*;x(t;y)) 

x(t) of system (1): ' — ^ ' ^^\, = . Let's construct the func- 

Pn+i{t;x{t;y)) Pn+i(0;y) 

tions 6i{t;x) = pi{t;x.) p~^^{t;x), I = l,n, taking into account the linear 

independence of the functions pi{0;y) and p„+i(0;y). 

Depending on the specific approach, described in the Introduction and 

the obtained independence of the functions 6i{t; x) with t = it was possible 

to construct a one-to-one correspondence: 

x, = qi{9i,e2,...,9n)- (34) 
Let's introduce the following notation just to shorten the entries: 

^ = {0i,02, ■ ■ ■ , 9n} , 

t = {qi,q2,---,qn} eW", 
^(t;x) = {pi(t;x),/32(t;x),...,p3(i;x)} GM*^. 
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Because of the condition (34) for some s > I and any other function xit] x) = 
x) p~^-|^(t;x) at the time moment t = can be represented hke this: 

x{0; y) = pn+s{0; y) Pnli (0; y) = ?(^(^)) • 

Then, consequently, for any function pn+s{t'-, x) being a kernel of the integral 
invariant (1) and because the lemma 1 conclusion, we get: 9i{t;x.): 

X(t;x(t;y)) = pn+sit;x{t;y)) p-l^{t;x{t;y)) = 

= pn+s{0; y) PnliiO; y) = x(0; y) = x(t; y)) • ^ (t; x(t; y)))^ . 

It follows from this relation the for any t > and each .s > 1: 

Pn+s{t;yi{t;y)) = p„+i(t;x(t;y)) V'(^'^(/3;;ji(i;x(t;y)) • (t; x(t; y))) ^ . 

Then the complete collection of the kernels consists of (n + 1) linearly inde- 
pendent functions. 

Corollary 1 The complete collection of linear independent stochastic first 
integrals of system (1) consists of n functions. 

Proof. The function ui{t;x) = pi{t;x)p~^-^^{t;x), possesses the proper- 
ties described in the definition 3 and presents the stochastic first integral of 
the system (1). There is n number of such functions. 

Remark 3 Let the random process x(t) is the solution of the generalized 
ltd 's SDE system that is presented like this: 

dtx{t) = a{t; x{t))dt + b{t; x{t))dw{t) + dP{t, A7) = dtx{t) + dP{t, A7). 

Then the generalized ltd 's formula (9) can be taken down as an equation: 

dtf{t; x{t)) = dtf{t- x{t)) + dtP{t, A7), (35) 

where dtf{t;x{t)) is the Ito's differential for the part dtx{t) and dtP{t, Aj) 
is the differential for the Poisson's component dP{t, A'j). 

Let's construct the equation for u{t; x) using the corollary 1 and the 
formula: 

\nus{t;x) = \nps{t;x) - In pi{t;x). (36) 
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By using (22), the Ito's generalized formula and (35) we'll do the differenti- 
ation of In p{t; x): 



dtlnp{t;x) 



1 



+ 



-dtp{t;x) 



1 



2p2(t;x) 



d{p{t;x)bik{t;x)) 



dt+ 



ln{ps (t;x - c/(t;x(t;y);7);7) • D (x ^(t;x;7)) } 
— lnps{t;x.) i'{dt;d'y). 



(37) 

Now let us examine the sum of the first two summands, dropping the func- 
tions' arguments: 



Si = - dtp - ^ 



d{pb, 



dxi 



dxi 



+ 



bjk) 



2 dxidxj 



1 



d{phi 



dt + 

2 

k) " 



dxi 



dwk{t) 



dxi 



dai 91np 1 d 

dxi * dxi 2p dx 
dbi k , dp^"^ 



dt 



dbikbjk , , , dp 
P r,^ + Oikbjk 



dxj 



dxi 



(38) 



2p2 y dxi dxi 
= S2dt 



, , , dlnp dbi k \ r / N 
dt - I bik— h -TT- I dwk{t) 



dxi 



dxi 



d\np dbik , . . 
Oik^T— + ^7— I dwkit). 



dxi 



dxi 
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We transform the part 

dai 



dxi * dxi dxj dxi 

2 dxidxj 2 ^ p dxidxj 
dbik 



dxi 



+ b 



dlnp\ [ dbjk , , dlnp\ 
i k o ) ( r Oj k 



daj _ 

dxj * dx 



1 d^jbikbjk) 

2 dxidxj 

I I dbik obi 

"2 



+ X bikbjk 



dxi J \ dxj dxj J 

dlnp ^ dbjkbjk dlnp ^ 

dxj dxi 
d^lnp 1^ ^ din p din p 

dxidxj 2 Qxj dxi 



dxi dxj 



+ 2bik'' ' 



dbj k dlnp din p din p 
~a fl ^ bikbjk 7^ — 

OXi OXi OXi OXi 



dxi 



dlnp dbikbjk dlnp 



dxi dxj dxi 
^ 1 d'\bikbjk) ^ 1 ^ ^ d'^lnp 1 dbjk dbjk 
2 dxidxj 2 *^ -^^dxidxj 2 SrE, 



dxj dxi 



(39) 



Hence, by substituting (39) into (38) and then into (37). We get the following 
equation for In p{t; x): 



dtlnp{t;x) = 
1 dbik{t;:x.) dbj kit; x 



dai{t;x) din pit; X 

- aiit;x) 



dxi 

+ ::bikit;x)bjkit;x) 



dxi 
d"^ lnpit;x) 



OXi dXj I oxiOXj 

dibik(t;x)bjkit;x)) din pit; x) 



+- 



dxi 



+ 



Id^ ibikit;x)bjkit;x)) 



2 

+ 



dxidxj 



bjkit;x) 



dxj 

dbikit;x) din pit; x) 



dxi 



dxj 



+ 



dt+ 



ln{ pit;x- git;xit;Y);-i);^) ■ D ix'-\t;x;-i)) } 
In pit; x) uidt; d'j) 



dbikit;x) dlnpit;x) 
-^-^+b,kit;x) 



dxj 



dwkit). 



(40) 



Relaying upon the obtained formula (40) we can construct the equations 
for psit;x) and piit;x), and having taken (36) into consideration we make a 
conclusion that the stochastic first integral uit; x; oj) of the generalized Ito's 
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equation is the solution for this equation: 
dtu{t; x) = 



-bikit;x] 



_d_ 

dxi 



, ,(?ii(i;x) 1 , , , .d'^u(t;:x.) 

— 7. h -bifc(t;x)6jfc(t;x) 

oxi 2 



6,fc(t;x) 



+ 



du{t; x) 
dxj 



dt 



u{t;x- g{t;x ^(t;x;7));7 



dxidxj 
u{t;x) v{dt]d'y). 



(41) 



The concept of the first integral for Ito's SDE system as the non random 
function over arbitrary random perturbation realizations has been introduce 
by the article [1]: 



ii(^0;x(0)j = u(^t;x(t;x(0)) 

We can also speak however about a stochastic first integral for Ito's SDE 
system. It should be observed that we can also speak about stochastic (but 
not only about determinate) first integral when Poisson's perturbations are 
missing. 

In case of just Wiener's perturbations (Ito classical equation) the in- 
vestigation of the properties of the first integral as a determined function 
was associated with stating the independence of such a function from re- 
alization w(i) []]. In the case under consideration from definition 3 the 

function M^t;x^ with its equation (41), depends over Poisson's perturba- 
tion are added (Ito generalized equation) such a requirement, as it follows 
from equation (41), leads to the following conditions. 

1 2 

Theorem 4 The random function u{t;:x.{t)) G C^'^ defined over the same 
probability space as the random process x(t) that is the solution for system 
(1) is the first integral of the system (1) iff the function u{t;x{t)) satisfies 
the terms C): 



1. bikit;x 



du{t; x) 



dxi 

perturbations ); 
^ du{t;x) ^ du{t;x) 



0, for all k = l,m (compensation of the Wiener's 



dt dxi 
of the time); 



ai{t;x) 



1, dbik(t;x) 
i^bjkit;x) 



dxi 



(independente 



3. u{t; x) — uyt;x + g{t;x;^)j =0 for any 7 G -^(7) in the whole field of 
the process definition (compensation of the Poisson's jumps). 
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Proof. Let's examine the equation (41), the sokition of which is the 

function u{t;x{t)). First let's carry everything to one side from the equals 

TT ■ 1 • 1 r / \ du(t; x) 

smg. Having taken mto account that dtu(t; xj = — — — at at we get for 



any t: 



du{t-x) du{t;x) 
+ ai(i;x)- 



dt 



+bikit;x] 



_d_ 

dxi 



bjkit;x[ 



dxi 
du{t; x) 

dxj 



dt 



1 , , , , , , 5^u(i;x) 
2&..(t;x)6,,(t;x)^-^+ 

dt + bik{t;x)^^^P^dwk{t)- 

OXi 



u(t;x- g{t;x ^{t;x;-f));j) + u{t;x) iy{dt;dj) = 0. 



R(7) 

Hence, the multipliers must be equal to zero with dt, dw{t) and v{dt]d'y). 
The Wiener's summand equaled zero: For the Wiener part we have: 

du{t; x) 



bik{t;:>i.)- 



dxi 



for all k = l,m. 



(42) 
(43) 



For Poisson's part we'll have: 

u(^t;x- g{t; x"^(t; x; 7)); 7^ - u{t; x) = 0. 

Now we re-arrange the equality (43) and proceed to variables: 

y = X - c/(t;x~^(t;x;7);7). 

But taking into account that x~^(t; x; 7) is the notation an inverse function 
for /(x) = X + (7(t;x;7), we get convinced that this term is equivalent to 
the following: 



n(t; x) — ti^t; X + (7(t; x; 7)^ = for any 7 G i?(7). 



(44) 



Then, using the rule of differentiating the product and the term (42), we 
obtain: 

d'^u{t;x) 



+ai{t; x) +-bik{t; x)bj k {t; x) - 



dxi 
d 



+- bik{t;x) 



d 



dx. 
du{t; x) 
dt 



, du{t; x) 
dxj 

du{t;x) 
+ ai(t;x) — h 



dxidxj 



, , .du(t;x) 



dxi 

9n(t;x) dbjk{t;x.) 



dx^ 



dxi 



du{t;x.) .du{t;x.) 1 du{t;x.) dbjk{t;x) 
^^ + "^(''")^^-2'^^^^*'")^^^ d^ 
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Consequently, 



du{t;x) du{t;x) r / sdbik{t;x 



dt ~^ dx 



ai{t;x) - - bjkit;x 



2 ' ^ dxj 



Thus, we get all the terms of C). 



Remark 4 // we analyzed the concrete realization then the non random 
function u{t; x) is the determinate first integral of the stochastic system. 

Remark 5 The concept of the stochastic first integral for a centered Pois- 
son measure was introduced in [3] . The obtained conditions for its realization 
take the necessity of determining the density of intensive Poisson distribu- 
tion. Thus, it makes no difference what is the probability distribution of 
intensities of Poisson jumps. This case is very important for constructing 
program controls [5]. 

The presented generalization of the Ito-Wentzell formula and the stochas- 
tic first integral concept [.'>] allow to construct of the program control with 
the probability is equaled to 1 for dynamic systems which being subjected 
to the Wiener perturbations and the Poisson jumps [">]. 
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